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Abstract. In this paper we consider spinor Bose-Einstein condensates with spin 
/ = 5 and / = 6 in the presence and absence of external magnetic field at the mean 
field level. We calculate all of so-called inert-states of these systems. Inert-states are 
very unique class of stationary states because they remain stationary while Hamiltonian 
parameters change. Their existence comes from Michel's theorem. For illustration of 
symmetry properties of the inert-states we use method that allows classification of the 
systems as a polyhedron with 2/ vertices proposed by R. Barnctt et at, Phys. Rev. 
Lett. 97, 180412 (2006). 
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1. Introduction 



Properties of Bose-Einstein condensates have been being found vastly since 1995 when 
first condensate was performed experimentallyp], |2l [31 IU [5] . After optical dipole traps 
were developed it was possible to create not only scalar but also spinor condensates. 
Computing properties of spinor condensates is more complicated than these of scalar one 
[U |5l [6]. In scalar condensates only one s-wave channel of interaction is possible since 
only one kind of bosons is present. On the other hand in spinor condensate of spin-/ 
atoms there are 2f + 1 internal states of atom. Therefore number of interaction channels 
increases. Because of this complexity finding stationary states for spinor condensates is 
not trivial even within the mean field approximation. 

In this paper we will first briefly show how to construct Hamiltonian for spinor 
condensates and how to perform calculations within mean field theory for that system 
[6]. After that we shall present Michel's theories and point out how it could help us 
finding stationary states [7]. The Michel's theorem ensures us that there have to exist 
so-called inert-states. Their name comes from a fact that they remain stationary for 
every values of Hamiltonian's parameters (such as external magnetic field, density etc.). 
Finally we will show how to use Michel's theorem and how to perform calculations 
of stationary states of spinor condensates [61 El |9]. All inert states for spin / < 4 
condensates and some inert states for higher spins have been presented in the literature 
0, [8j [9]. We will present all inert states for spin-5 and spin-6 condensates with and 
without external magnetic field. For illustration of symmetry properties of the states 
we use spherical harmonics expansion and the method proposed by Barnett et. al 



2. Hamiltonian and mean field theory 



That part is very well described in [6] so we shall briefly point way for constructing 
Hamiltonian and its most important assumptions. In spinor condensates atoms possess 
internal degrees of freedom characterized by a state of hyperfine structure. We 
should remember that usually in spinor Bose-Einstein condensates description by / we 
understand total angular momentum of all electrons and nuclei of atom which we will 
just call spin. So the wave function of spinor condensate in the mean field description 
should be written as: 



N 



=1 



(1) 



*-/(r<) 

where \I/ m (rj) is wave function component describing particles with projection m of spin 
on chosen quantization axis. 

Let us denote ^ m {r) as a bosonic field operator which creates a boson at point 
r with spin projection equal m. In the absence of external trapping potential but in 
the presence of external magnetic field, energy of this system can be decomposed into 



Inert-states of spin-5 and spin-6 Bose-Einstein condensates 3 

kinetic energy, energy of Zeeman effects (linear and quadratic) and interaction energy. 
So we can write: 

H = H + V , (2) 

where by H we understand part of Hamiltonian without interaction between particles. 
Assuming the direction of external magnetic field B along the z-axis it can be 
decomposed as follows: 



Ho = J d 3 f £ ^ 



ft 2 v 2 

■—^-pUz) m\m 2 



c'2\ 



2M 



^m 2 {r) , (3) 

in 1 ,m 2 =-f 

where M is mass of atoms, p ~ B factor describes the linear Zeeman effect, q ~ B 2 
factor describes the quadratic Zeeman effect and (f z ) mi m 2 is spi n matrix (in the chosen 
coordinates system it is diagonal (f z ) mi m 2 — m i^m 1> m 2 )- 

In dilute and ultra-cold atomic gases two-body s-wave interaction channels are 
important only. From wave function's symmetry considerations we can find out that the 
only possible channels for identical bosons correspond to even values of total spin F of 
two interacting particles. So the interaction part of the Hamiltonian can be decomposed 
in the following way 

V= £ . (4) 

F=0,2,...,2/ 

We can construct (using Clebsh-Gordon coefficients) operators which annihilate pairs 
of bosons with total angular momentum F and its projection M on the z-axis 

/ 

A FM (r,r')= £ (F,M|/,m i; /,m 2 )^ mi (f)^ m2 (f') . (5) 

m 1 ,m 2 =-f 

Now we can see that operators can be written as: 

\M = I fd'r f dW F \f,f') £ A^ M (f,f')A FM (r,r') . (6) 

* J J M=-F 

Function (r, r ') describes spatial dependence of interactions between two particles 
in the channel F. This function using s-wave scattering length a F can be assumed to 
be: 

v W{r,r') = g F 8®{r-r') , (7) 
where g F = ^jj-a F . 

We can also use two relations to simplify our final form of the Hamiltonian. Namely: 
• Hilbert space completeness relation: 

F 

E E A ] FM {r)A FM {r)=:n{r)n{r): , (8) 

F=0,2,...,2f M=-F 

where : : stands for normal ordering of the creation and annihilation operators, 
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relation coming from composition of angular momentum 

F 



/(/ 



A FM {?)A FM {r) 



:F(r)-F(r):= £ £ 

F=0,...,2/ M=-F 

In the particle number conserving version of the mean field theory we have to 
calculate mean value of the Hamiltonian assuming that N bosons are in the same state, 
see Eq. ([TJ. We assume also that the spinor and spatial parts of the wave- function 
factorize, i.e. \^{r)) = </>o(r)|0. In the absence of an external trapping potential the 
ground state mode 0o = where V is the volume of the system, and the energy 
reduces to the expression 



E(C) = {(\H\() 



where \() 



7^ (E f m =-f (m°ln;o) 1°) ■ Vector C = (Cf, -;C-F) 
parameter. This parameter completely describes state of spinor condensate. Neglecting 
terms proportional to 1/N final result for the energy per particle reads [6] 



N 



(10) 

) T will be called an order 



N 



f 

E 



-/ 



ICr, 



+ E 

F=0,...,2/-4 



—pm + qm 
d F p j2 



:Pc 



;pci\F\ 



Afm-Afm 



M=-F 

where p is the particle density. Parameters cq, c\, di are linear combinations of 
coefficients gi and are presented in Appendix 1. Functions used in (ITT]) have the following 
forms 

/ 



mi ,rri2 = 



-/ 



.4 



FM 



J 

E 

mi ,7712= 



(F,M\f, mi;/, m 2 )CmiCn 



(12) 



(13) 



In order to calculate all stationary states of the energy ( ITT]) one has to derive the 
corresponding Euler-Lagrange equations and solve them which usually has to be done 
numerically. However, a certain class of stationary states can be found analytically by 
consideration of symmetries of the system. These are the so-called inert-states which 
we present in the following sections. 



3. Michel's theorem [7] 

Michel's theorem allows us to find special stationary states. They are called inert-states 
because they are stationary for any values of Hamiltonian's parameters (in our case the 
parameters are: particle density, scattering lengths and external magnetic field). The 
stationarity of these states is guaranteed by symmetries of a system. 

When external magnetic field is present the symmetry group of the energy function 
( ITTj) is a product group Gb^o = U{1)f z x U(l). Part U{1)f z is rotation about the field 
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direction (i.e. z-axis) and part U(l) is just global phase changing symmetry. If an 
external magnetic field is not present the symmetry group is Gb=o = 5*0(3) x U(l). 
Meaning of the part U(l) is the same, i.e. it is the global phase changing symmetry, 
whereas SO (3) is group of rotations in 3-dimensional space. 

For our further considerations we should call symmetry group by G and by e(Q 
function defined in a finite-dimensional smooth domain Ai which has this symmetry. 
By orbit G(Q of point ( we understand set of points being results of action of all group 
elements on point ( 

G(0 = {gC\g eG}<zM. 

Isotropy group of point £ is a set of group elements which act on Q as an identity element 

G< = {ge G\g( = (} c G . 
We say that any two isotropy groups are conjugate if 

3geG:G ( = g-G c -g- 1 . 
Stratum is union of all orbits of points whose isotropy groups are conjugate 

S(0 = IJ G(C) c M . 

We say that an orbit is isolated in its stratum if there exists neighborhood of this orbit 
U D G(() so that its intersection with stratum gives that orbit 

G(() = u n S(C) . 

Michel proved following theorem. 

Theorem 1 (Michel) Every function defined in a finite- dimensional smooth domain 
A4 and having real numbers as co- domain which is invariant under elements of compact 
Lie group G i.e.: 

ye : M -> 3? : V# G G A V( G M : e(C) = e(g() , 

have common orbits of extrema which are orbits being isolated in their stratum. 

From this theorem we know that for our spinor condensate system we are able to 
find some stationary points ( just by finding isolated orbits. The crucial point of the 
Michel's theorem proof is to show that gradient of function e(Q is a vector which is 
tangential to the stratum. We know from definition of orbit that gradient of function is 
orthogonal to the orbit. And if any orbit is isolated in its stratum we know that gradient 
has to be null on it so this orbit has to be extremal (minimal, maximal or saddle) point. 

Barnett et. al. [TO] showed how to translate symmetries of spin-/ states \() into 
symmetries of a polyhedron with If vertices. We should not redo it here but point 
out the most important steps. We have to find a set of maximally polarized states, i.e. 
f\x) = f\x) j which are orthogonal to \(), 

(C|X>=0. (14) 
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The states \x) can be parameterized as x = e^tan | where and 9 denote position of 
vertices of a polyhedron on the unit sphere. For our considerations the most important 
corollary is that maximal number of found vertices is equal to 2f because Eq. (I14j) is 
actually a polynomial equation of the 2f degree. Therefore when we are looking for 
inert-states of the system we may consider only these point group whose corresponding 
regular polyhedron has no more than 2f vortices. 

Now we should present way of calculating and finding inert-states [61 13 19]. In order 
to do that we should go through all subgroups of a global symmetry group of a system. 
For each subgroup we should find states which are invariant under this subgroup. In 
that case this subgroup is their isotropy group. 

First we will show how to calculate inert-states for subgroups with only one 
generator. In our case such generators have form K = e 1 * ^ (i.e. they are rotation 
operator about f2-axis by an angle 0). It is worth pointing out that if K belongs to the 
symmetry group of the system so is e lX K because in our case the symmetry group of the 
system is a direct product of a rotation group and the global phase changing symmetry 
group. All eigenvalues of K can be written as e n . If we find all eigenvectors of a 
generator K we can identify some of them to be inert-states. Through eigensubspaces 
of K we should choose these which are non-degenerate. That means that for each 
eigenvalue e n there exists only one eigenvector £ 7 . Therefore we know that group 
generated by e~ %1 K is isotropy group of this vector. Moreover we see that there is only 
one orbit whose isotropy group is this group. Thus, this orbit must be isolated in its 
stratum because stratum consists of this orbit only. Then, the Michel theorem tells us 
£ 7 is an inert-state. 

If some subgroup has two or more generators, as in the case of point groups, we 
ought to find vectors which are common eigenvectors for all generators. If there is only 
one common eigenvector it is an inert-state of the system. When on the other hand 
the common eigenspace is not one-dimensional we still can calculate stationary states 
which are not inert-states. To this end we have to find extrema of the energy of the 
system within the common eigensubspace. Existence of such stationary states comes 
from another Michel's theorem j6]. However, these states depend on the Hamiltonian 
parameters. We shall not consider them in this paper. 

4. Inert-states for spin-5 condensates with external magnetic field 

Global symmetry group of spinor condensates with external magnetic field is direct 
product group Gb^o = U(l)p z x U(l). General group element has form: 

<7(0,7) = e*e-^, (15) 

where <fi is global phase change and 7 is angle of rotation about the i?-axis which coincides 
with the direction of the magnetic field. We have found 11 inert-states for this case. 
All of them are inert-states whose isotropy groups are continuous. Analysis of discrete 
subgroups does not reveal any additional inert-state. 
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Figure 1. Visualization of the state (f5+, Eq. (TIT)) , with Barnett et al. method [TU] . 
Red circles indicate positions of 2/ = 10 vertices of a polyhedron. In the presented 
case all vertices are situated on the same point (that is indicated by the number 10 in 
the figure) and the polyhedron is reduced to a single point. 




k 2tz 

Figure 2. Visualization of the state Cf5+, Eq. (jTT|) . by means of spherical harmonics. 
Surface \^(9, 4>)\ 2 = 1, see Eq. (fT6|) . is plotted in the figure. Colours express phase of 
the state at a given point. 

We present the inert-states with their isotropy group's generator and the 
corresponding energy. Isotropy groups are denoted by U(l)F z + a (<f>)- Only for one state 
we have decided to illustrate its symmetry. One way is that proposed by Barnett et. al. 
|10j . Second one is by using spherical harmonic expansion. That is, every bosonic spin 
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state can be expand in terms of spherical harmonics as follows 

9(0, 4)= £ Cm>7(M) (16) 

m=-f 

We can visualize states by painting surfaces |\1/(6 I ,0)| 2 = 1. Colour of surface should 
express phase of state at a given point [101 E2] • 

• U{l) Fz+ ^ 
generator: {e^e - * 9 ^} 

Cf 5+ = (1,0, 0,0, 0,0, 0,0, 0,0, Of (17) 

,. , 1 25 
^(Cfs+) = ^pc + yPC! + 25q - 5p , (18) 

The state is visualized in Fig. [T] and Fig. [21 

• U(l) Fz+4<t> 
generator: {e Al ^e~ l ^ z } 

( F4+ = (0,1, 0,0, 0,0, 0,0, 0,0, Of (19) 
£(Cf4+) = \pc Q + 8p Cl + 16q - Ap , (20) 

generator: {e 3 ^ - "^} 

C ra+ = (0,0, 1,0, 0,0, 0,0, 0,0, Of (21) 

1 9 14 

£(Cfs+) = ^pc + ^pci + —pd 6 + 9q - 3p , (22) 

• U(l) Fz+2<t> 
generator: {e 2l ^e _l ^ z } 

Cf 2+ = (0,0, 0,1, 0,0, 0,0, 0,0, Of (23) 
^(Cf2+) = \pc Q + 2p Cl + J^pd 4 + ^pd 6 + 4q-2p, (24) 

• U(l) Fz+(P 
generator: {e^e - *^} 

C^i + = (0,0, 0,0, 1,0, 0,0, 0,0, Of (25) 

. * 1 1 25 , 10 , 14 , 

e(Cfi+) = + 2 pCl + 286 P 2 + 143 9 4 + 187 P 6 + q ~ P ' ^ ^ 



£/(lk 

generator: {e~*^ z } 



Cp = (0,0,0,0,0,1,0,0,0,0,0) 7 (27) 

. , 1 1 , 25 , 9 , 40 , 

= + ^ + ^pd 2 + ^4 + ggjpde , (28) 
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generator: {e _t ^e~ l ^ z } 

Cfi- = (0,0, 0,0, 0,0, 1,0, 0,0, Of (29) 

, . 1 1 25 , 10 , 14 , 

e(Cfi-) = ^pc + -pci + 7^P«2 + Y^ pdi + I87 prf6 + g + p ' ( 30 ) 

• ^(1)^-20 

generator: {e~ 2lrt> e~ lrt>fz } 



Cf2- = (0,0, 0,0, 0,0, 0,1, 0,0, Of (31) 

1 35 , 84 

2 PC0 + 2/?Cl + 286^ 4 + 935' 



1 35 84 

e(( F 2-) = -pc + 2pci + — pd 4 + w^pde + M + 2p , (32) 



^(1)^-30 

generator: {e _3 ^e~*^ z } 



C ra - = (0,0, 0,0, 0,0, 0,0, 1,0, Of (33) 

1 9 14 

2 PC ° + 2 pCl + 85' 



1 9 14 
e(Cra-) = tjP c o + oP c i + ^P^6 + 9g + 3p , (34) 



^(l)F z -40 

generator: {e _4 ^e _l ^ z } 



Cf4- = (0,0, 0,0, 0,0, 0,0, 0,1, Of (35) 
s(Cf4-) = l^pco + 8pci + 16g + 4p , (36) 

generator: {e -5 *^ - * 9 ^} 

Cfs- = (0,0, 0,0, 0,0, 0,0, 0,0, If (37) 
1 25 

£(Cfs-) = -pco + ypc 2 + 25g + hp . (38) 

5. Inert-states for spin-5 condensates without external magnetic field 

In this case global symmetry group is Gb=o = 5*0(3) x U(l). General group element 
can be written as: 

g{a, (3, 7 , 0) = e^e^'e'^e^' , (39) 

where is global phase change and a, (3, 7 are Euler angles. This group has the same 
continuous subgroups as the global symmetry group considered in the previous section. 
So we expect that the corresponding inert-states are the same. It easy to show that in 
the case without magnetic field energies of some of these inert-states are equal. That 
is, for p = q = we get e(( F5+ ) = e(( F 5-), ^(Cf4+) = e(Cf4-), ^(Cf3+) = £(Cf3-), 
s((f2+) = £((f2-) and e((fi+) = ^(Cfi-)- We can show that members of each pair 
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can be transformed to each others by applying some elements of the group Gb=o- So 
they belong to the same orbit and therefore they are actually the same states. We have 
decided not to show these states one more time but present how to change one vector 
into another by action with group elements: 

(7(0,7t,0,0)Cf5 + = Cf5- (40) 

^(0,7T,0,7t)Cf4+ = CF4- (41) 
3(0,7r,0,0)CF3 + = CF3- (42) 
^(0,7r,0,7r)CF2 + = CF2- (43) 

3(0,7t,0,0)Cfi + = Cfi- • (44) 

Now we should consider all discrete point groups. Let us begin with the dihedral 
group D n . Generators of this group are C njX (rotation about z-axis by angle 2n/n) 
and C 2)X (rotation about x-axis by angle ti with simultaneous phase change by ir). We 
would like to stress that in our consideration these generators will be supplemented with 
additional phase factors. This is necessary in order to obtain proper isotropy groups and 
we are allowed to do that because the symmetry group of our system contains global 
phase changing freedom. It is worth mentioning that the dihedral group with infinite n 
is also an isotropy group of the previously found state (p, Eq. ( 127)) . In the limit n — > oo 
the dihedral group becomes continuous. Other inert-states for the dihedral group are 
presented here: 

generators: {e l7r C 2 , x , e l7T C WjZ } 

Cmo = 0, 0, 0, 0, 0, 0, 0, 0, 0, -1) T . (45) 

We can show that C^zCdio = e~ ln (oio and Cio^Cdio = g~ Z7T (dw- Thus, in order 
to obtain the isotropy group we have to supplement the original dihedral group 
generators with the phase factor e l7T . 

. , 1 1 , 75 , 9 , 15 , 

s(Cd W ) = ~ 2 Pc, + ^Pd + W2 pd, + — pdt + —pd 6 . (46) 

generators: {e t7r C 2 ,x, e t7r C StZ } 

Cmo = ^=(0, 1, 0, 0, 0, 0, 0, 0, 0, -1, 0) T (47) 

1 1 3 9 96 

e(Cr)io) = + 7^P d o + J^P^ + 143^ + 935^ ' ^ 48 ^ 
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D 6 

generators: {e t7r C 2tX , e t7r C^ z } 

1 



( m = -==(0, 0, 1, 0, 0, 0, 0, 0, -1, 0, Of (49) 

111 9 2689 

e(Om) = ~ 2 Pc, + ^Pd + ^Pd 2 + —pd A + —pde , (50) 



generators: {e t7r C 2tX , e^C^} 

1 



Cm = ^(0, 0, 0, 1, 0, 0, 0, -1, 0, 0, Of (51) 

1 1 3 9 96 

e(Cm) = -pc + — pd + -^Pd2 + j^Pd* + ~^P de ' ( 52 ) 

• D 2 

generators: {e l7T C 2 , x , e 47r C 2 , z } 

Cd2 = ^(0, 0, 0, 0, 1, 0, -1, 0, 0, 0, Of (53) 

111 9 41 

e(Cm) = ~pc + — pd + —pd 2 + j^pd* + ^P rf 6 • (54) 

Another point group, which should be considered is tetrahedron group T with 
generators C 3)X+y+z (rotation about x + y + z-axis by angle 2n/3) and C 2;Z (rotation 
about 2- axis by angle tt). There is only one inert-state whose isotropy group is this 
group, namely: 

Ct = ^(0, 1,0, 2, 0,0, 0,-2, 0,-1, Of . (55) 
One can check that 

C 3>x+y+z ( T = e i2 r( T , (56) 
C 2 , Z ( T = Ct • (57) 

Thus the generators of the isotropy group are {C 2jZ , e~ l27r ^C 3>x+y+z }. 

/ j. \ 1 6 , 15 , 168 , 

s(Ct) = - 2 PC + — Pd 2 + — pd 4 + —pd 6 . (58) 

The state with tetrahedral symmetry is visualized in Fig. [3] and Fig. [4j 

There are not any states with the octahedral group as isotropy groups. We know 
also that maximal number of vertices of a polyhedron which represents states for spin-5 
condensates is 10. Since icosahedron has 12 vortices there cannot be any inert-state 
having the icosahedral group as its isotropy group. Therefore the states presented are 
all inert-states one can find for spin-5 condensates. 
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z 




Figure 4. The same as in Fig. [2]but for the state (t, Eq. (|55|) . 

6. Inert-state for spin-6 condensate with external magnetic field. 

Here we have the same global symmetry group as for the corresponding spin-5 case. 
The calculations are similar as before and results are the following: 

generator: {e^e - *^} 

£ F6+ = (1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, Of (59) 

£(£f 6+ ) = ipc + 18pci + 36q - 6p , (60) 

• U(l) Fz+5(j} 
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generator: {e 5i<?i e _ ^ z } 

& 5+ = (0,1, 0,0, 0,0, 0,0, 0,0, 0,0, Of (61) 
x 1 25 

e(fF5+) = 2 pc ° + Y pCl + 25q ~ 5p ' ^ 



U{1)f x +a4> 

generator: {e 4i *e _i ^ z } 



£f4+ = (0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, Of (63) 

1 45 

-OCn + 8pci H / 

T p 1 266' 



1 45 
s{£ F4+ ) = -pc + 8pc! + — pd 8 + 16q - Ap , (64) 



generator: {e 3i ^e _i ^} 



£ F3+ = (0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, Of (65) 
1 9 42 12 



generator: {e 2 ^e~*^ z } 



= (0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, Of (67) 

. 1 245 252 , 18 , 

£{£f2+) = -pc + 2pci + ^^P«4 + ggggMs + 2~47 P 8 + 4 9 - 2p , (68) 

U(1)f m+4 > 

generator: {e^e - *^} 

Zfi+ = (0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, Of (69) 

. , 1 1 21 , 140 , 210 , 180 , 
e(M = 2^+2 PCl+ ^ 2+ 243T pd4+ 3^ p4+ 2W pd8+g ^ ' (70) 

U(1)f. 

generator: {e"^ 2 } 



^p = (0,0, 0,0, 0,0, 1,0, 0,0, 0,0, Of (71) 

, , 1 1 7 , 126 , 200 , 175 , , , 

£(£p) = ~ 2 PCo + Y6 pd + W3 pd 2 + — P d A + —pd 6 + —pd 8 , (72) 



generator: {e _t ^e~ l ^ z } 



= (0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, Of (73) 

, . * 1 1 21 140 , 210 , 180 , 

s{Cfx-) = 2 PC0+ 2 PCl+ ^ 2+ 243T prf4+ 3^ P4+ 2W Prf8+g+P ' (?4) 
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generator: {e -2 *^ - * 9 ^} 

£ F2+ = (0,0, 0,0, 0,0, 0,0, 1,0, 0,0, of 

, . \ 1 n 245 , 252 , 18 , A n 
e{^F2-) = 7>pc + 2pci + T^yP^ + ^^P rf 6 + 7^jP d s + ^q + 2p 

generator: {e -3 *^ - * 9 ^} 

£ F3+ = (0,0, 0,0, 0,0, 0,0, 0,1, 0,0, of 

1 9 42 12 

e(ff3-) = + 2 pCl + 323 prff3 + I33 pd8 + 9q + 3p ' 



generator: {e~^e~ l ^^ z } 



T 



£ F4 - = (0,0,0,0,0,0,0,0,0,0,1,0,0) 

1 45 

e(f j?4-) = oP c o + 8p c i + T^P^s + 16g + 4p 
z zoo 



generator: {e -5 "^ - *^*} 



£ F5 _ = (0,0, 0,0, 0,0, 0,0, 0,0, 0,1,0) 

1 25 
e(fF5-) = + yP c i + 25 <? + 5 P 



generator: {e _6 ^e~*^ z } 



^ 6+ = (0,0,0,0,0,0,0,0,0,0,0,0,1 
e(£ F6+ ) = ^pc + 18pci + 36g + 6p . 

7. Inert-states for spin-6 condensates without external magnetic field. 

Similarly as in the corresponding case of spin-5 condensates the inert-states found in 
the presence of an external magnetic field are also inert-states without the field. Also 
similarly as previously pairs of these states are actually the same states because in the 
absence of the field they can be transformed to each other by applying elements of the 
global symmetry group of the system. These transformations are showed below: 

^(O, tt,0, 0)fre+ = fre- (85) 
(?(0,7r,0,7r)£ F5+ = £ F5 _ (86) 



<?(0,7r,0,0)^ 4+ = £43- (87) 
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3(0,7r,0,7r)O, 3+ = £F3- (88) 

0(O,7r,O,O)fo+ = fr2- (89) 

(7(0,7r,0,vr)e F1+ = eFi- • (90) 

The state £p, Eq. ( ITT]) , is also an inert-state whose isotropy group is the continuous 
dihedral group D^. The other inert-states with the dihedral groups as isotropy groups 
are the following: 

generators: {e l7r C 2 ,x, e 47r Ci 2 , z } 

W = 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1) T (91) 

. , 1 1 , 11 , 891 , 11 , 11 , . , 
e(M = 2^ + ^ + ^ 2+ 1^^+646^+6916^ ' < 92 > 

generators: {e I7r C 2 ,x, e I7r Ci , 2 } 

£mo = i(0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, Of (93) 

1 1 11 99 275 275 

e(foio) = 2^ + ^ 0+ ^ 2 + ^ 4+ 2584^ + 6916^ 8 ' (94) 

generators: {e t7r C 2 , x , e t7r C 8 , 2 } 

£ D8 = i(0, 0, 1, 0, 0, 0, 0, 0, 0, 0, -1, 0, Of (95) 

, 1 1 , 11 , 1152 , 8 , 3589 , 
s(U) = 2^+26^0+ 1^+17^+3^53^+ 19019^ ' ( 96 ) 

generators: {e l7r C 2iX , e J7r C 6i2 } 

£ D6 = -L(0, 0, 0, 1, 0, 0, 0, 0, 0, -1, 0, 0, 0) T (97) 

t , N 1 1 , 25 , 729 , 3697 , 3793 , , 
e(tm) = 2^+2^^°+4^^ 2 +3lb^^ 4 +2^^ 6 +7^^ 8 ' (98 

• D 4 

generators: {e l7r C 2>a; , e l7T C^ z } 

{Ida = ^=(0, 0, 0, 0, 1, 0, 0, 0, -1, 0, 0, 0, Of (99) 
, , 1 1 , 25 , 537 , 373 , 601 , , , 

= 2 p ^ + Y6 pdo+ Tm pd2+ !^ pd4+ ^ pde+ ^ pd8 ' (100) 
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D 2 

generators: {e^C^, e %n &2,g} 

t D2 = -L(0, 0, 0, 0, 0, 1, 0, -1, 0, 0, 0, 0, Of (101) 

/ \ 1 1 , 99 , 1002 , 155 , 655 . . 
e{tm) = 2^ + ^^ + Tb^^ 2+ T^^ + ^^ + 19W9^ ^ 

In the case of the tetrahedron group we have found two candidates for inert-states: 
generators: {C 2 , 2 , e- i27T/3 C 3:X+y+z } 



,y/TT i y/E iy/7 y/E i y/TT. T 

f T1 = (1 — 0, — =, 0, — , 0, 0, — , 0, — =, 0, - — T , 103 
8 Ay/2 8 4 ' ' 8 Ay/2 8 J v ; 

/ j. \ 1 72 , 15 , 224 , 2922 , . 

e{M = - 2 Pc» + ^pd, + — pd 4 + — pd 6 + —pd 8 , (104) 

• generators: {C 2 ,z, e l27T/3 C 3 , x+y+z } 

£t 2 = (-g-, 0, - i - 7 =, 0,— ,0, -_, 0, -, 0, - i7 =, 0, _) , (105) 

x 1 72 , 15 224 , 2922 , . . 

e{M = -pc + — pd 2 + — pd 4 + — prf 6 + — P d 8 . (106) 

The energies of these two states are the same and it turns out that they belong to the 
same orbit 

3(0,0,7r,7r)£ T1 = fr 2 , (107) 

and thus they are actually the same inert-state. 

Another group which have to consider is the octahedron group. The octahedron 
group is group generated by three generators C& z (rotation about z-axis by angle 7r/2), 
C 3)X+y+z (rotation about x + y + z-axis by angle 27r/3) and finally C 2yX+y (rotation about 
x + y-axis by angle it). For this group we have found one inert-state: 

to = (0, 0, 0, 0, 0, 0, 0, 0, 0, 0) T , (108) 

, 1 1 , 1323 4 , 189 , , , 

e(to) = ~ 2 pc + -pd + — prf 4 + — pd» + —pd 8 . (109) 

The to state is an eigenstate of the C*4 jZ , C^+j^ and C 2 , x + y generators corresponding 
to eigenvalues equal 1. Thus, no phase supplement is necessary because the original 
generators form the isotropy group of to- 

The last group is the icosahedron group. We ought to consider that group because 
an icosahedron has 12 vertices and maximal number of vertices for illustration of spin-6 
states is exactly 12. Icosahedron group's generators are (% )2 (rotation about i?-axis by 
angle 2ir/5), C 2)V (rotation about y-axis by angle it) and C 2< (-i+^5) x +( 2 +^E) z (rotation 
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about (-1 + y/5)x + (2 + y/E) i?-axis by angle n). For this group we have found one 
inert-state: 

£y = ^(0, - V7, 0, 0, 0, 0, 0, 0, 0, 0, y/7, 0) T , (110) 

1 1 121 . . 

= -pco + ^gP d + —p d 6 . (Ill) 

Similarly as in the case of the octahedron group no phase supplement is necessary 
because the generators C^ z , C^y and C , 2 ,(-i+ v / 5)x+(2+ v / 5)z f° rm t ne isotropy group of £y. 



8. Conclusion 



In this paper we have presented all inert-states for both spin-5 and spin-6 systems with 
and without external magnetic field. We have found all these states by considerations 
of symmetries of energy density functional. Idea of finding inert-states was used in 
superconductors and Bose-Einstein condensates for lower spins before [131 EH El El E] • 
Some inert-states for spin-5 and spin-6 were also presented in [9]. 

We look for the inert-states employing the Michel's theorem. Yip argues that if 
one has found inert-states for some group he need not to consider any subgroup of this 
group because there are no additional inert-states related to the subgroups [8]. It is 
obviously truth. We want to pay attention, that any of dihedral groups presented in 
our paper is subgroup of other. However, Makela and Suominen in Ref. [9] point out 
that we should not forget about very important global phase changing factors in their 
generators, which cause that they cannot be subgroups of other. They consider groups 
_D 4 and D 8 . These groups without phase factors obviously obey relation D 8 C D 4 . But 
in our case it is not true, because: 

(e^C 8 , z ) 2 = C 4 , z ^ e w C M . 

Thus we can easily see that we cannot construct elements of group D± by using group 
D 8 elements if they are supplemented by phase factors [9] . Consequently both _D 4 and 
D 8 has to be considered separately in investigation of inert-states of the systems. The 
same argumentation can be repeated for groups D 2 and _D 4 . 

The results presented in this paper can be used as a first step for calculating all 
stationary states of considered systems. Changing Hamiltonians parameters (external 
magnetic field or scattering lengths) one can also identify quantum phase transitions 
when ground states of the systems change their symmetries. Dynamical properties of 
the phase transitions can be analyzed employing the Bogoliubov approach that allows 
us to investigate stability properties of the stationary states [T5l HEl [17]. In order to 
analyze stability of the inert-states one has to solve the Bogoliubov-de Gennes equations. 
These equations, however, can not be solved analytically and numerical calculations are 
necessary. Because we do not know values of scattering lengths for any spin-5 and spin-6 
condensates numerical studies of the stability of the inert-states are postponed for future 
research. 
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Appendix 1 

Values of coefficients of Hamiltonian for spin-5 system. 



Co = 


^(25# 8 - 6g w ) 




Cl = 






d = 


Y^( 19 9o + 36£fio 


- 550 8 ) 


d 2 = 


^(1902 + 330x0 


- 520 8 ) 


di = 


^(1904 + 26^0 


- 450 8 ) 


d 6 = 


^(1906 + 150x0 


- 340 8 ) 



Values of coefficients of Hamiltonian for spin-6 system. 



Co = 


1 . 

7^(360io - 130x2) 




Cx = 


^(012 — 0io) 




d = 


1 , 

— (230o + 550l2- 


780io) 


d 2 = 


^(2302+520 1 2- 


750i O ) 


G?4 = 


^(2304+450l2- 


680 1O ) 


da = 


^(2306 + 340x2 - 


570io) 


d 8 = 


^(230 8 + 190i2- 


420io) 
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